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ABSTRACT 

This  paper  presents  results  of  an  analysis 
of  ground  shock  data  from  more  than  one  hundred 
explosion  tests  conducted  In  soil  over  the  past 
35  years.  Burst  positions  varied  from  fully 
buried  to  contact  detonations  in  soil  and  for 
shallow  depths  Into  concrete  protective  overlays. 
Soil  conditions  ranged  from  loose  dry  sand  to 
saturated  clay.  Empirical  equations  are  pre- 
sented that  predict  the  magnitude  and  time  his- 
tories of  the  expected  stresses  and  ground 
motions  as  a function  of  burst  position,  soil 
indices  and  burster  layer  thickness , 

BACKGROUND  \ 

More  than  one  hundred  explosion  teats  have 
been  pertormed  over  the  past  35  years  In  soil  to 
characterize  the  ground  shock  produced  by  bulled 
conventional  munitions.  Lampson  (1)  conducted 
pioneering  tests  during  World  War  II  with  burled 
charges  in  several  soil  types  from  which  cube 
root  scaling  was  developed  and  verified.  Follow- 
ing Lampoon's  lead,  the  MOLE  and  UET  tests  (2,3) 
were  conducted  in  the  early  1950' s to  extend  the 
ground  shock  data  base  to  other  soils  and  rocks. 
Since  that  time  many  other  burled  burst  tests 
have  been  conducted  in  a wide  variety  of  soils 
with  charges  from  1 pound  to  S'-O  tons  (4-11), 

Many  other  small-scale  teats  have  goue  unreported. 

Testa  performed  for  special  projects  pro- 
duced the  first  measurements  of  stresses  directly 
beneath  near-ourface  explosions  on  burster  slabs 
(12,13).  Coupling  of  shallow  explosions  from 
penetrating  weapons  was  investigated  by  Ingram  (8) 
in  the  CENSE  program  and  by  Drake  and  Little  (un~ 
reported).  Brown,  et  al  (14),  investigated  the 
propagation  of  ground  shock  through  rock  rubble 
overlays.  All,  totaled,  more  than  50  tests  have 
been  conducted  just  to  determine  explosive  cou- 
pling from  partially  buried  weapons  in  burster 
slabs . 

Ground  shock  data  from  these  tcBts  were 
dram:  together  and  analyzed  to  provide  an  update 
to  the  Army  design  manual,  TM  5-855-1,  "Funda- 
mentals of  Protective  Design  for  Conventional 
Weapons,"  The  result  of  that  analysis  is  the 
basis  for  this  paper.  Hue  to  the  limited  space 
provided,  only  the  empirical  prediction  equa- 
tions are  given. 


GROUND  SHOCK  THREAT 

The  ground  shock  produced  by  bombs  exploding 
on  or  within  the  ground  near  buried  structures  gen- 
erally provides  the  dominant  threat  to  these  facil- 
ities. Stresses  from  buried  bursts  can  be  greater 
in  magnitude  and  much  longer  duration  than  corre- 
sponding bursts  in  air. 

Significant  enhancement  of  the  stresses  and 
ground  motions  will  occur  as  the  weapon  penetrates 
more  deeply  into  the  surrounding  Boll  or  backfill 
before  it  denotates.  Often,  protective  layers  of 
concrete  or  rock  rubble  are  provided  over  the 
structure  to  limit  the  weapon  penetration  thus  re- 
ducing the  effective  coupling  of  the  explosion  and 
Increasing  the  weapon  standoff. 

Important  variables  affecting  the  intensity  of 
the  loading  are:  a)  weapon  size  and  distance  to 

the  structure,  b)  the  mechanical  properties  of  the 
soil  or  rock,  and  c)  the  de^t’u  of  penetration  of 
the  weapon.  Of  these  factors,  the  effect  of  soil 
or  rock  properties  is  least  predictable  by  simple 
handbook  methods.  Ground  shock  intensity  may  vary 
by  more  than  two  orders  of  magnitude  when  the  soil 
is  varied  from  low  density  dry  soils  to  saturated 
clays . 

There  are  two  important  cases  to  consider  in 
assessing  the  ground  shock  threat  to  burled  facili- 
ties: 1)  bombs  that  explode  overhead,  generally 

on  or  within  the  protective  concrete  or  rock  rub- 
ble overlay,  causing  a direct  loading  of  the  roof 
slab,  and  2)  weapons  that  penetrate  into  the  sur- 
rounding soil  and  detonate  beside  the  facility 
loading  the  walls  and  floor.  While  the  same  gen- 
eral ground  shock  prediction  equations  apply  for 
both  cases,  the  role  of  the  site  geology  and  the 
protective  overlay  requireo  a somewhat  different 
application  of  these  equations.  These  cases  are 
shown  in  Figure  1 . 

SOIL  1R0PERTY  EFFECTS 

Ground  shock  propagation  in  earth  media  is  a 
complex  function  of  the  dynamic  constitutive  prop- 
erties of  the  soil,  the  explosive  products  and 
the  geometry  of  the  explosion.  No  single  soil 
index  or  combination  of  Indices  can  adequately 
Jescribe  thi.B  process  in  a simple  way  for  all 
cases. 

Water  can  tiave  a profound  influence  on  ground 
shock  propagation  in  cohesive  soils,  particularly 
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Figure  1.  Geometry  for  explosion  against  a 
buried  facility 

as  saturation  reaches  95  percent  or  greater. 
Typically  water  is  fully  bound  with  the  skeletal 
structures  for  these  soils,  providing  signifi- 
cant contribution  to  the  overall  stiffness  and 
strength  of  the  soil  structure.  As  saturation 
approachs  100  percent , pronounced  increases  in 
peak  stresses  and  accelerations  have  been  ob- 
served in  wet  clays,  clay  shales  and  sandy  clays. 
Stresses  similar  to  shock,  waves  in  free  water 
have  been  noted  In  saturated  clays.  Saturation 
as  measured  by  free  standing  water  in  boreholes 
may  not  be  an  accurate  measure  of  the  true  sat- 
uration depth,  particularly  where  seasonal  water 
table  fluctuations  introduce  small  amounts  of 
air  Into  the  soil.  Seismic  surveys  generally 
will  show  a sharp  Jump  in  the  wave  speed  to  more 
than  5000  feet  per  second  at  this  depth. 

Granular  soils  with  high  relative  density 
are  generally  not  strongly  influenced  by  water 
saturation  as  are  cohesive  soils.  The  stiffness 
f granular  soil  is  provided  by  the  grain  to 
grain  contacts  in  the  skeleton  with  only  a small 
contribution  by  the  free  water.  Consequently, 
controlled  laboratory  and  field  experiments  in 
dense  nearly  saturated  sands  did  not  show  large 
influence  of  the  pore  water  on  the  resulting 
shock  wave  propagation.  However,  the  effects  of 
W'-ter  in  low  relative  density  sands  can  produce 
similar  effects  as  those  seen  in  cohesive 
soils  [15].  In  these  cases,  the  soil  skeleton 
can  collapse,  and  the  grain  to  grain  contact 
lost  resulting  in  high  pore  pressures  as  the 
sand  liquifies.  These  sites  would  not  normally 
be  considered  for  construction  of  hardened 
facilities . 

Seismic  velocity,  c , is  often  used  as  a 
crude  index  of  noli  or  rock  properties  for  ground 
shock  prediction  purposes.  It  provides  a simple 
measure  of  the  stiffness  and  the  density  of  the 
soli  thru  the  relationship 


where  M Is  the  stiffness  or  modulus  of  the  soil 
and  p is  Its  mass  density.  The  seismic  veloc- 
ity also  provides  a relationship  between  distance 
and  time. 


A great  dual  of  caution  must  be  used  in  gen- 
eralizing the  use  of  seismic  velocity  as  a ground 
shock  index.  Cementation  in  granular  soils  such 
as  dry  desert  alluvium  may  result  in  abnormally 
high  propagation  velocities  (approaching  4000  feet 
per  second).  Yet  these  materials  may  exhibit  very 
high  air  filled  voids  and  low  relative  densities — 
qualities  that  would  classify  them  as  very  poor 
transmitters  of  ground  shock.  Low  seismic  veloc- 
ities, on  the  otner  hand,  would  generally  indicate 
poor  ground  shock  transmission  qualities. 

The  attenuation  rate  with  range  of  the  ground 
shock  magnitude  is  controlled  by  the  irreversible 
crushing  of  the  void  volume  within  the  soil  matrix 
by  the  passage  of  the  stress  wave.  In  cohesive 
soils,  the  volume  of  the  air  filled  voids  is  an 
index  for  aLtenuatlon  of  ground  motions,  while  the 
beat  index  for  attenuation  rates  in  granular  media 
is  the  relative  density.  Because  relative  density 
is  not  always  available,  dry  unit  weight  can  be  an 
effective  index  for  ground  shock  attenuation. 

Soils  with  high  relative  density  (high  dry  density) 
or  low  air  voids  will  attenuate  the  ground  shock 
more  slowly  than  low  relative  density  or  nigh  air 
void  materials.  Figure  2 shows  the  peak  stress 
from  explosions  in  typical  soils. 


Figure  2.  Peak  stress  from  explosions 
in  varying  soil  types 


GROUND  SHOCK  PREDICTIONS 

Stress  and  particle  velocity  pulses  can  be 
cha-ructerized  by  exponential  like  time  histories 
that  decay  rapidly  in  amplitude  and  broaden  as 
they  propagate  outward  from  the  explosion.  The 
characteristic  time  for  these  time  histories  can 
be  measured  in  arrival  time  from  the  source, 
t . where 


R is  the  distance  from  the  explosion  and  c i9 
the  seismic  or  propagation  velocity.  Typically 
these  wave  forms  rise  sharply  to  the  peak  with 
the  riae  time,  t 
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this  is  about  1/10  of  the  travel  time  to  the  tar- 
get point.  From  the  peak.,  the  pulse  decays  mono- 
tonically  with  lime  to  nearly  zero  in  1 to  3 
travel  times  given  by  the  following  equations 


P(t)  - F e 


V(t)  - V fl  - Bt/t  )e 
o'  a 


t >_  0 (3a) 


t > 0 (3b) 


where  p(t)  is  the  stress,  V is  the  particle 
velocity,  and  a and  8 are  time  constants. 
While  the  time  constants  generally  vary  with  spe- 
cific site  conditions,  they  can  be  taken  as 


for  moat  applications,  PQ  and  VQ  are  values  of 
the  peak  stress  and  perLicle  velocity  to  be  deter- 
mined by  the  following  equations.  Other  waveform 
parameters  such  as  impulse,  displacement  and  ac- 
celerations may  be  derived  from  these  functions. 

Since  the  characteristic  time  is  Inversely 
proportional  to  seismic  velocity,  explosions  in 
high  velocity  media  such  as  saturated  clay  will 
produce  very  short,  high  frequency  pulses  with 
high  accelerations  and  low  displacements.  On  the 
other  hand,  detonations  in  dry  louse  materials 
will  cause  much  longer  duration,  low  frequency 
ground  motions. 

Peak  particle  velocity  and  peak  stress  are 
related  by 


where  p is  the  mass  density.  Free-fleld 
stresses  and  ground  motions  from  bombs  detonating 
on  and  within  burster  layers  or  In  the  soil  along 
side  the  structure  ate  given  in  the  following 
expressions : 


where  ro  J.s  the  peak  pressure  (psi) , V0  is  the 
peak  particle  velocity  (ft/sec),  yQ  is  the  peak 
acceleration  (g's),  dQ  is  the  peak  displace- 
ment (ft),  I0  is  the  peak  Impulse  (psi-sec),  f() 
Is  the  mass  density  (ib-sec/ f t^) , c is  the  seis- 
mic velocity  (ft/sec),  pc  is  the  acoustic  imped- 
ance (psi/ft/sec) , R is  the  distance  to  the  ex- 
plosion (ft),  W is  the  charge  weight  (lb),  and  f 
is  the  coupling  factor  for  near  surface  detona- 
tions. For  preliminary  design  considerations  the 
following  table  is  suggested  for  selecting  the 
seismic,  velocity,  acoustic  Impedance  and  attenua- 
tion coefficients: 

SUCCESTED  COEFFICIENTS  FOR  DESIGN 


Seismic 

Velocity 

c 

Acoustic 

Impedance 

U'C) 

Attentuntion 
Coef  f icient 

Material  Description 

ft/sec 

psi/ft/sec 

n 

Loose,  dry  sands 
and  gravels  with 
low  relative  density 

600 

12 

3-3.25 

Sandy  loam,  loess, 
dry  3ands  and 
backfill 

1000 

22 

2.75 

Dense  sand,  with 
high  relative 
density 

1600 

44 

2.5 

Wet  sandy  clay  with 
air  voids  (greater 
th**n  h percent) 

1800 

66 

2.5 

Saturated  sandy 
clays  and  sands 
with  small  amount 
of  air  voids  (less 
th3n  1 percent) 

5000 

48 

2.25-2.5 

Heavy  saturated 
clays  and  clay 
shales 

>5000 

150-180 

1.5 

A more  detailed  description  is  provided  in 
Table  1 for  soils  encountered  in  explosion  test  pro- 
grams. Simple  soil  parameters  such  as  wet  and  dry 
unit  weights,  air  filled  voids  and  seismic  velocity 
are  shown  to  ansist  in  relating  the  explosion 
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effects  parameters  to  the  design  soil  conditions. 
Note  that  the  attenuation  coefficient  and  seis- 
mic velocity  are  closely  related  to  dry  unit 
weight  for  granular  soils  while  air  void  con- 
tent is  important  for  cohesive  soils. 

GROUND  SHOCK  COUPLING  FACTOR 

The  magnitude  of  the  stress  and  ground  mo- 
tions will  be  greatly  enhanced  as  the  weapon 
penetrates  more  deeply  into  the  soil  or  the  pro- 
tective burster  layer  before  it  detonates.  A 
concept  of  an  equivalent  effect  coupling  factor 
is  introduced  to  account  for  this  effect  on  the 
ground  shock  parameters  and  is  defined  as  follows: 

The  coupling  factor,  f , is  defined  as 
the  ratio  of  the  ground  shock  magnitude 
from  partially  to  shallow  buried  weapon 
to  the  ground  shock  magnitude  from  a 
fully  buried  burst  in  the  3ame  medium. 

f (P»  V»  d,  I,  a)  near  surface 
(P,  V,  d,  I,  a)  contained 

A single  coupling  factor  is  applicable  for  all 
ground  shock  parameters  that  depends  upon  the 
depth  of  burial  of  the  center  of  the  weapon  and 
the  medium  being  penetrated,  i.e.,  soil,  concrete 
or  air.  It  is  important  to  note  that  the  coupling 


factor  concept  used  her  a does  no?;  produce  an  equiv- 
alent charge  but  rather,  "it  is  scale  factor  to 
reduce  the  ground  shock  computer  iron  a burled 
burst  with  the  full  charge  weight  .0  account  for 
the  shallow  burial. 

Coupling  factors  are  different  for  bursts 
in  air,  soil  and  concrete  and  depend  upon  the 
scaled  depth  of  burst  of  the  weapon.  These  fac- 
tors are  shown  in  Figure  3.  The  coupling  factor 
for  air  is  a constant 


0,14 


and  is  recommended  for  contact  bursts. 

In  the  case  where  a weapon  penetrates  into 
more  thaw  one  material,  i.e.,  a long  bomb  that 
penetrates  the  concrete  slab  and  is  partly  buried 
in  soil,  the  coupling  factor  is  computed  as  the 
sum  of  the  coupling  factors  in  each  of  the  mate- 
rials weighted  in  proportion  to  the  charge  weight 
contained  within  each  medium. 


■ 1 &) 


(6) 


where  f is  the  total  coupling  factor,  is  the 

coupling  factor  for  each  component  material,  i.e., 
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air,  soil,  concrete,  Is  the  weight  of  the 

charge  In  contact  with  each  component  material, 
and  W is  the  total  charge  weight.  Since  most 
bombs  are  cylindrical,  the  coupling  can  also  be 
defined  as 

f ■ 2 fi (r)  <5 * 7 

where  is  the  length  of  the  weapon  in  contact 

with  each  material  and  L is  the  total  weapon 
length. 

DISCUSSION  AND  CONCLUSIONS 
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Empirical  expressions  were  derived  from  a 
fit  to  a large  body  of  ground  shock  data  from 
buried  and  near-surface  bursts  in  soil.  Several 
Important  observations  were  made  concerning  the 
role  of  soil  properties  on  scaling  of  ground 
shock: 

1.  Near  the  explosive  source,  peak  parti- 
cle velocities  in  soils  «*end  to  a single  curve 
thet  is  nearly  independent  of  the  soil  proper- 
ties. This  observation  can  be  explained  in  part 
by  the  interaction  of  the  detonation  wave  in  the 
explosive  with  the  soil. 

2.  Peak  stresses  scale  In  proportion  to 
the  seismic  velocity. 

3.  Attenuation  of  the  peak  ground  shock 
magnitudes  is  strongly  dependent  on  the  rela- 
tive density  in  granular  soils  or  to  the  air 
void  volume  In  cohesive  soils.  Because  the  seis- 
mic velocity  is  also  Influenced  by  these  param- 
eters, the  attenuation  coefficient,  n , can  be 
estimated  from  the  seismic  velocity  as 
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c(fps) 


500-600 

750-1000 

1000-1400 

1400-1800 

>5000 


3-3.5 

3 

2.75 

2.5 

1.5-2.25 


4.  Time  scales  in  proportion  to  the  time  of 
arrival.  Thus,  the  pulse  tends  to  spread  in  pro- 
portion to  the  distance  traveled,  with  a rise 
time  of  about  1/10  of  the  travel  time  and  a dura- 
tion on  the  order  of  2-3  travel  times. 

5,  Because  of  the  time  scaling,  peak  accel- 
erations are  proportional  to  the  seieraic  veloc- 
ity, peak  displacement b are  inversely  proportional 
to  the  seismic  velocity  while  the.  peak  impulse  is 
only  sensitive  to  density  variation. 
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TABLE  1 SOIL  PROPERTIES  FROM  EXPLOSION  TESTS 


Soil  Description 

Dry  desert  alluvium  and  playa, 
partially  cemented 

Loose,  dry,  poorly  graded  sand 

Loose  wet  poorly  graded  sand- 
free  standing  water 

Dense  dry  sand,  poorly  graded 

Dense  wet  sand,  poorly  graded- 
free  standing  water 

Very  dense  dry  sand,  relative 
density  =1002 

Silty-clay,  wet 

Moist  loess,  clayey  sand 

Wet  sandy  clay,  above  water- 
table 

'Saturated*  sarrd-below  water 
table  in  marsh 

'Saturated'  sandy  clay,  below 
water  table 

'Saturated1  sandy  clay,  below 
water  table 

Saturated  stiff  clay  saturated 
clay-shale 


(1)  High  because  of  cementation. 

(2)  Estimated. 


i 

!.i 


Dry  Unit 
Weight 

Ydry 

lb/ft3 

Total  Unit 
Weight 

Air-Filled 

Voids 

% 

Seismic 

Velocity 

c 

ft/sec 

Acoustic 

Impedance 

pc 

psi/f t/sec 

Attenuation 

Coefficient 

n 

87 

93-100 

>25 

2100-4200  ^ 

40 

3-3.25 

80 

90 

>30 

600 

11.6 

3-3.5 

97 

116 

10 

500-600 

12.5-15 

3 

99 

104 

32 

900-1300 

25 

2.5-2.75 

108 

124 

9 

1000 

22 

2.75 

105 

109 

30 

1600 

44 

2.5 

95-100 

120-125 

9 

700-900 

18-25 

2.75-3 

100 

122 

5-10 

1000 

28 

2.75-3 

95 

120-125 

4 

1800 

48 

2.5 

— 

— 

1-41'23 

4900 

125 

2.25-2.5 

78-100 

110-124 

1-2 

5000-6000 

130 

2-2.5* 

100 

125 

<1 

5000-6600 

130-130 

1.5 

..M. 

120-130 

0 

>5000 

135 

1.5 

